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To study the electron-phonon mechanism of superconductivity in Na^OoC^ • 1/H2O, we perform 
semiquantitative analysis of the electron-phonon interaction (EPI) between relevant optical phonons 
(breathing and shear phonons) and t^ a electrons (a\ g and e' g electrons) in the presence of trigonal 
distortion. We consider two kinds of contributions to the EPI; the EPI originating from the Coulomb 
potential of O ions and that originating from the d-p transfer integral between Co and O in C0O6 
octahedron. We find that the EPI for shear phonons, which induces the interorbital hopping of 
electrons, is large in Na^CoC^ • J/H2O because of the trigonal distortion of C0O2 layer. For this 
reason, T c for s-wave pairing is prominently enlarged owing to interorbital hopping of Cooper pairs 
induced by shear phonons, even if the top of e' g electron band is close to but below the Fermi level 
as suggested experimentally. This mechanism of superconductivity is referred to as the valence- 
band Suhl-Kondo (SK) mechanism. Since the SK mechanism is seldom damaged by the Coulomb 
repulsion, s-wave superconductivity is realized irrespective of large Coulomb interaction U ~ 5 eV 
at Co sites. We also study the oxygen isotope effect on T c , and find that it becomes very small due 
to strong Coulomb interaction. Finally, we discuss the possible mechanism of anisotropic s-wave 
superconducting state in Naa;Co02 • J/H2O, resulting from the coexistence of strong EPI and the 
antiferromagnetic fluctuations. 
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I. INTRODUCTION 

Superconducting layered cobalt oxide Na^CoC^ ■ 2/H2O 
(x ~ 0.35, y ~ 1.3) with T c ~ 4.5 K has attracted consid- 
erable attention since it is the first discovered supercon- 
ducting cobaltates [l[. The parent compound Na x Co02 
reveals rich phase diagram depending on the Na content 
x 0. Yokoi et al. have found that the electronic state 
of Naa;Co02 drastically changes at the boundary x ~ 0.6 
3. For example, the uniform magnetic susceptibility de- 
creases with decreasing temperature for x < 0.6, while 
it shows a Curie- Weiss-like behavior for x > 0.6 [1, 0|- 
The quasiparticle spectra in photoemission spectroscopy 
Q and optical conductivity in infrared spectroscopy Q 
also decrease with decreasing temperature for x < 0.6. 
It is noteworthy that these "weak pseudogap" behavior 
is also observed in the normal state of superconduct- 
ing Na^OoC^ • 2/H 2 0. In the superconducting state of 
Naa;Co02 • 2/H2O, a sizable decrease of the Knight shift is 
observed independently of the direction of magnetic field, 
which indicates that the spin singlet superconductivity is 
realized ■ A power law behavior in \ jT{F suggests 

the large anisotro py in superconducting gap, like a line- 
node state [Tol. ITll.Tl2| . On the other hand, the decreasing 
rate of T c due to non-magnetic impurities in this system 
is much smaller than that of the d-wave superconductor, 
and it is as small as that of s-wave superconductor MgB 2 
[il^ . This fact suggests that the sign of the superconduct- 
ing gap function is unchanged everywhere. 

In considering the mechanism and the pairing sym- 
metry of superconductor, the Fermi surface (FS) topol- 



ogy gives the most important information. The FSs in 
Na^CoC^ are composed of £2 9 -orbitals of Co ion, which 
split into the ai s -orbital and twofold e^-orbitals due to 
crystalline electric field. The first principle calculation 
based on local density approximation (LDA) [l4| had pre- 
dicted the presence of six small hole pockets due to the e' g 
bands near the K-points, in addition to a cylindrical FS 
due to the a\ g band around T-point. However, such e' g 
hole pockets are not observed in ARPES measurements, 
since the e' g bands are completely below the Fermi level 
independently of the Na content [l5[ . They are also not 
observed in the bulk-sensitive ARPES using soft X-ray 
that has longer escape depth Moreover, the shape 
of FS docs not change by the intercalation of water [Tj] • 
We stress that an estimated specific heat coefficient 7 es t 
using the a\ g Fermi velocity given by ARPES [3 is con- 
sistent with experimental vaule 7 ~ 11 mJ/mol K 2 . If 
there were e' g hole pockets, the realized 7 should be about 
3 times greater than the experimental value because of 
the large density of states (DOS) given by the e' g hole 
pockets, as pointed out in Ref. [Hj]. This fact reinforces 
the observation by ARPES measurements. 

To determine the FS topology theoretically, present 
authors have studied the normal electronic state of 
Naa;Co02 based on the fluctuation-exchange (FLEX) ap- 
proximation, which is a self-consistent spin-fluctuation 
theory (l9| . In the FLEX approximation, experimentally 
observed weak pseudogap behaviors in the DOS and mag- 
netic susceptibility appear when the top of the e^-band is 
slightly below the Fermi level, that is, the e' g hole pockets 
are absent. We found that the weak pseudogap bchav- 
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iors originate from (i) antiferromagnctic (AF) fluctua- 
tions due to Coulomb interaction, and (ii) large DOS of 
the e' g hole pockets that exist slightly below the Fermi 
level. When e' g hole pockets are present, on the other 
hand, ferromagnetic fluctuations are induced by the e' g 
hole pockets, and " anti-pseudogap behavior" appears in 
the DOS These results are highly inconsistent with 
experiments. Therefore, Na :c Co02 should have a single 
cylindrical FS around the T-point. 

After the discovery of Na.2,Co02 • 1/H2O, various kinds 
of superconducting states had been proposed 0, [2l], ■ 
In particular, possibility of triplet superconducting state 
due to Coulomb interaction were investigated based on 
the FS with e' g hole pockets, by using perturbational 
theory [H, 13] or FLEX approximation [25|, How- 
ever, successive experimental efforts had revealed that e' g 
hole pockets are absent [TtJ , and singlet superconducting 
state is realized 0, Moreover, ferromagnetic fluctu- 
ations are not observed by inelastic neutron diffraction 
in a superconducting single crystal Therefore, we 

now have to find a mechanism of singlet superconducting 
based on the large a\ g FS around the T point. If e' g hole 
pockets are absent, however, the expected magnetic fluc- 
tuations seem to be too small to realize unconventional 
superconductivity (25[. Considering the small impurity 
effect on T c [lU, we should examine the possibility of 
superconductivity caused by the electron-phonon inter- 
action (EPI). In fact, several experimental studies show 
the presence of considerable electron-boson coupling in 
Naa;Co02: For example, the kink structure was observed 
in the quasiparticle spectrum in ARPES measurement at 
uj ~ 600 cm -1 [l6j]. In addition, w-dependence of scatter- 
ing rate in infrared spectroscopy shows a steep upturn 
at u = 500 ~ 600 cm -1 [6j. These experiments suggest 
the strong electron-boson coupling. Since the energy of 
these bosonic modes corresponds to the frequency of rel- 
evant optical phonons, EPI is expected to be strong in 
Na x Co0 2 . 

We have studied the electron-phonon mechanism for 
superconductivity in our previous paper [28j . by notic- 
ing two kinds of the optical modes (breathing (Ai g ) and 
shear {E g ) modes) that strongly couple with 3c? electrons 
in Co. We have found that the shear mode phonon, 
which represents the oscillation of O ions parallel to the 
C0O2 layer, induces the transition of Cooper pairs be- 
tween different bands. Due to this mechanism, a con- 
siderably strong pairing interaction for s-wave pairing is 
realized, which is known as the Suhl-Kondo (SK) mech- 
anism [23]. The SK mechanism works even if e' g hole 
pockets are absent, as long as the top of the e^-band 
(valence band) is close to the Fermi level compared with 
phonon frequency [28j |. In Ref. f28j . we have discussed 
that s-wave superconductivity due to this "valence band 
SK effect" is expected to be realized in water-intercalated 
Naa;Co02 • yH 2 since the top of the valence band is sup- 
posed to approach the Fermi level [1(3, [HI . We note that 
extended s-wave scenario had be en prop osed based on 
the two concentric a\ g FSs model |32|, [33J . 



However, we did not take account of the large Coulomb 
repulsive interaction in 3c? orbitals of Co in Ref. [28[, 
which works to prevent the s-wave pairing. It is an im- 
portant problem to elucidate whether or not the attrac- 
tive force due to EPI overcomes the Coulomb repulsion 
and an s-wave superconductivity can be realized. For 
that purpose, we have to know the values of EPI and 
the Coulomb interaction precisely. According to a re- 
cent first-principle cluster calculation based on a quan- 
tum chemical ab-initio method, the Coulomb interaction 
at Co site in Na x Co02 is 4 ~ 5 eV [13] ■ Also, the mass 
enhancement factor in M2,Co02 (M=Na, K or Rb) due to 
optical phonons at ~ 0.1 eV is estimated to be m* /m ~ 2 
by ARPES measurement [TB]. However, there is little in- 
formation about the precise matrix elements of EPI in 
Na a ,Co0 2 . 

In this paper, we quantitatively examine the EPI be- 
tween 3c? electrons and relevant optical phonons that in- 
volve the deformation of CoOg octahedron. The EPI 
originates from both the change of the Coulomb poten- 
tial and that of the transfer integral due to the displace- 
ment of O ions. We calculate the value of EPI using the 
second order perturbation of the transfer integral. We 
find that the EPI for shear phonon is considerably in- 
creased by the trigonal distortion of C0O2 layer, which 
had not been taken into account in our previous study 
[28| . Thus, s-wave superconductivity can be realized in 
Naa;Co02 • J/H20 irrespective of the large Coulomb in- 
teraction U ~ 5 eV. We also study the oxygen isotope 
effect on T c , and find that it becomes very small when 
U = 4 ~ 6 eV. 

This paper is organized as follows. In ClTI we derive 
the EPI microscopically in case with and without distor- 
tion of C0O2 layers. We discuss the change of EPI due 
to intercalation of water. The erroneous result of EPI 
given in Ref. [28| is corrected. In mill we explain the 
3-band model for ti g electron system, and we derive the 
linearized gap equation. In ijlVI numerical results of the 
transition temperature are obtained by solving the gap 
equation numerically. In SjVl we discuss the robustness 
of s-wave superconductivity over the strong Coulomb in- 
teraction. We also discuss the oxygen isotope effect and 
the possible mechanism of anisotropic s-wave supercon- 
ducting state in Naa;Co02 ■ 2/H2O. Finally, results of this 
paper are summarized in flVIl 



II. ELECTRON PHONON INTERACTION IN 

Na.CoOj yH 2 

In this section, we derive the EPI in Na x Co02 ■ 2/H2O 
microscopically. This compound consists of two dimen- 
sional C0O2 layers separated by a thick insulating layer 
composed of Na ions and H2O molecules. C0O2 layer 
comprises a triangular network of CoOg octahedra that 
share edges as shown in Fig. [1] (a). In this network, Co 
ions (black circles in Fig. H](a)) and O ions in the upper- 
layer (white circles) and lowerlayer (gray circles) forms 
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triangular lattices, respectively. Due to the crystalline 
electric field by octahedral coordination of O ion, five- 
fold 3d orbitals of Co ion split into threefold ti g (d xy , d yz 
and d zx ) orbitals and twofold e g (d x 2_ y 2 and d 3z 2_ r 2) 
orbitals. Hereafter, we neglect e g -orbitals since they are 
completely empty. In Na^OoO^, C0O6 octahedra are 
trigonally distorted along c-axis, so that C0O2 layers be- 
come thinner. Due to this trigonal distortion, threefold 
ti g orbitals split into the a\ g orbitals and twofold e' g or- 
bitals as shown in Fig. [2l According to the LDA calcula- 
tions, the a\g orbital forms large hole-like FS around the 
T-point and the e' g orbitals form six hole pockets near 
the K-points. Since Na ions and O ions are monova- 
lent cation and divalent anion, respectively, valence of 
Co ions is Co^ 4 ~ x ^ + (mixed valence of Co 3+ and Co 4+ ). 
Therefore, the number of 3d electron is 5 + a; ,that is, the 
number of hole in ti g orbitals is 1 — x in Na^Co02 ■ 2/H2O. 
If the concentration of oxonium ion (H30 + ) is z, the hole 
number is given by 1 — x — z [35[ . Since the present pa- 
per is devoted to the s-wave superconductivity due to the 
EPI, obtained results depends not on the hole number, 
but on the top of the e^-bands measured from the Fermi 
level. For this reason, we do not consider the existence 
of the oxonium ion hereafter. 

Next, we determine the optical phonon modes that 
cause strong EPI. In Na^CoOo, there are 14 kinds of 
optical modes at zone center [361 ] . Since the potentials of 
tig orbitals are considerably changed by the deformation 
of CoC>6 octahedron and each C0O2 layer is almost iso- 
lated, we have only to consider single Co02-laycr phonon 
modes. According to Ref. [28[, there are 4 single-layer 
phonon modes. Among them, the ungerade modes in 
which all O ions move in the same directions has no cou- 
plings with tig electrons within the linear term of the 
displacement. As a result, only the A\ g mode (breathing 
mode) and the E\ g mode (shear mode) strongly couple 
with tig electrons. In breathing mode and shear mode, O 
ions oscillate in the direction parallel and perpendicular 
to the C0O2 layer, respectively, as shown in Fig. Q] (b). 

Here, we describe the method of calculation of the EPI. 
We focus on the zone center modes and calculate the EPI 
via the frozen phonon method. To obtain the EPI, we cal- 
culate the change of potentials of ti g orbitals due to the 
displacement of six O ions in CoOg octahedron, which 
is shown in Fig. [T] (c). We consider two types of con- 
tributions to the EPI; one is the Coulomb potential by 
2 ~ ions, and the other is the effective potential due to 
the transfer integrals between Co and O. In Ref. [28[, 
we consider the EPI due to transfer integrals up to the 
fourth order processes. However, the third order term 
was wrong in sign. After correcting this mistake, we ver- 
ified that a reliable EPI can be obtained by the second 
order perturbation since the third order term and the 
fourth order term almost cancel. In the present study, we 
find that the effect of trigonal distortion on EPI, which 
was not taken into account in the previous study p8l |. is 
significant. The calculation of the EPI (in the presence 
of trigonal distortion) will be shown in Al fi jll Bl) in 



detail. 




FIG. 1: (a), (b) Crystal structure of a C0O2 layer. Co and O 
in upside (Oi, O2, O3) and downside (O4, Os, Oe) of C0O2 
layer forms triangular lattice, respectively, (c) The structure 
of C0O6 octahedron without trigonal deformation. Note that 
z-axis does not correspond to c-axis. 




trigonal 
distortion 

FIG. 2: Crystalline electronic field for d-electron in Naa;Co02. 
The splitting between a\ g and e' g is caused by the trigonal 
deformation of CoOe octahedron. 



A. EPI in the absence of trigonal distortion 

Here, we calculate the EPI between ti g electrons and 
optical phonons in the absence of trigonal distortion of 
C0O6 octahedron. We set the xyz coordinate system 
as shown in Fig. [1] (c). c and a crystal axes are along 
the (1,1,1) and (1,-1,0) directions in the xyz coordinate 
system, respectively. The coordinates of six O ions are 
(±a,0, 0), (0, ±a,0), (0,0, ±a), where positive (negative) 
sign corresponds to Oi, O2 and O3 (O4, O5 and Oe) ions 
that locate upper (lower) side of C0O2 layer. In breath- 
ing and shear modes, O ions oscillate in the direction 
parallel and perpendicular to the c crystal axis, respec- 
tively. The displacement due to a breathing phonon is 
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u br _ _|_tL_^, 1,1), where positive (negative) sign cor- 
responds to the displacements of Oi, O2 and O3 (O4, O5 
and 06 ) ions. The displacement due to shear phonon 
is expressed as u SH = ±(vi,V2,— v\ — v%). Here, we 
choose the orthogonal bases as u SH1 = ±^^=-(1, — 1, 0) 

and u SH2 = ±^=-(1,1,-2), where positive (negative) 
sign corresponds to the displacements of Oi, O2 and O3 
(O4, O5 and 06 ) ions. 
Here, we calculate the EPI originates from the 



Coulomb potentials by using the point charge model. 
The Coulomb potentials Vc(r) of an electron with charge 
— e at the center of six O 2- ions is expressed as Vc(r) = 

Ei=i,6^1F^7T' where r * =1 " 6 ) arc thc P ositi on 
vectors of O ions. Owing to the displacement of O 
ions due to a-mode phonon, Vc(r) changes to V§{r) = 

Si=i 6 | r _ ri 1 _ u a| ■ To calculate thc EPI, we expand 
8Vq(y) = Vq(t) — Vc(r) up to the first order in u a as 
follows: 



8Vg*(r) = I + ^(xy + yz + zx)) u™, (1) 

47T£o \ a or I 

sv ° m v = H (-^ - y>) - ~ ">) " s "' (2 » 

Wc ' H2(r) = + v '- 2 ' 2) + i? [2xy - y '- zx} ) " SH2> (3) 



where we have dropped thc fourth and higher order terms 
in x, y and z. Thc matrix elements of thc EPI due to 
Coulomb potential is given by (^\5Vq\i>) , where fx, v = 
x,y,z. Since the wave function for /ii/-orbital (fj,v = 
xy,yz, or zx) is expressed as (f)^ v {r) = vTo"^r0(r), it 
is easy to show that (yz\xy\zx) = \{r 2 ), where (r 2 ) 



I (j> 2 (r)d 3 r is the expectation value of the square of the 
radius in 3d-orbital. We use the notation r 2 — (r 2 ) here- 
after. Similarly, (zx\yz\xy) = (xy\zx\yz) = jr 2 ,. Thus, 
we obtain the EPI for breathing phonon as 



SV 



BR 



2e 2 

47T£ 



1 
1 
1 



4V3r 2 



i i 

- - 

7 y 7 

- - 

7 7 u 



,BR 



where the first, the second and the third column (row) 
correspond to d xy , d yz and d zx orbitals, respectively. 

Here, we change the basis of t2g-orbitals from the (d xy , 
dyz, d zx )-hasis into the (a\ g , e'g 1 , eg 2 )-basis, which is the 
basis of the irreducible representation under the trigonal 
crystalline electric field. The transformation matrix U 
from the (d xy , d yz , c? Z2 )-basis to the {a\ g , e'g 1 , eg 2 )-basis is 
given by 



\ai g ) 
K) 

\J2\ 




(4) 



By operating U and W on the left side and the right side 
of the matrix in Eq. ([?]), respectively, the matrix form of 



SV 





e' 1 e' 2 
c g > c g 


BR 




C = 


(! 






2e 2 


2V3 


47re 


a 2 


2e 2 


2V3 


47re 


a 2 





a? 
a 2 c 



u BR , 



2Mu B r 



2A/WBR 



(5) 



1 _ 7af ' 
1 + gf)- 00 



Here, wc introduce a dimcnsionlcss parameter u a — 
u a I y/h/2Mu> a that represents a displacement due to 
a-modc phonon as thc unit of zero-point motion 
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y/H/2Mu> a , where ui a is the frequency of a-mode phonon. 
Then, u a is expressed by the sum of creation and anni- 
hilation operators for the a- mode phonon, u a =b ol + b\ l . 
Next, we calculate the EPI for shear phonons. Using 



(xy\x 2 \xy) = (zx\x 2 \zx) = fr^ and (yz\x 2 \yz) = ±r 2 , 
we obtain the matrix form of the EPI for shear phonons 
from Eqs. 10 -© in the (d xy , d yz , d zx )-basis as 



SV 



SHI 



SV 



SH2 




3V2r 2 



SHI 



V6r 2 



SH2 



SHI 



(8) 



,SH2 



(9) 



By transforming Eqs 

'9 ' C S 



(dig, e'g, e' 2 )-basis, we obtain 



into the expression in the 



SHI 



SV 



~ SH2 

sv c = 






b? 










-b'i 


- 


-b$ 











-b? 





"2 





-b? 





-b§ 



u sm , 



u sm , 



(10) 
(11) 



b? 



(«i s l^c SH1 | 
2e 2 2y/3 



47T£o I 
„/l| <n^SH2i 



Km 



c 



2Mw SH 7a 2 ' 
; 1 )=6f/(4V2). 



(12) 
(13) 



Thus, the EPI originating from the Coulomb potential 
is represented by four coupling constants af , , frf and 
6?- 

Next, we derive the EPI that originates from the 
changes of d-p transfer integrals between Co and O. For 
this purpose, we consider the CoOg octahedron in which 
Co 4+ ion (hole number is unity) is surrounded by six 
2 ~ (no holes) ions. In this octahedron, the uncertainly 
principle allows the virtual process in which a hole in d 
orbitals transfers to p orbitals and turns back. By this 
second-order process, the effective potential of a hole in 
d orbitals is lowered by — t 2 d /A p d, where t p d is the d-p 
transfer integral and A p d(> 0) is the charge transfer en- 
ergy. When the transfer integral is changed by St p d, the 



effective potential of a hole becomes ~(t p d + St p d) 2 / A p d- 
Thus, in the electron representation, the change of the 

effective potential is given by SVr — 2t p dSt p d/ A p d- In 

y a 



2 »}' 





Oi 

s- — ■ 



a 



'pdn 



(X) 



FIG. 3: Transfer integrals between d xy orbital in Co and 
(p x ,Py) orbitals in O. In this figure, d xy electron can transfer 
to p y orbitals in Oi and O4 via ±t p d7r. 



C0O6 octahedron, the total change of the effective po- 
tential (8Vt)w for d-orbital is given by 
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(5V, 



T )H' 



E E 

i= 1 — 6 m—p x ,p y ,p z 



t pd (i, I, m)5t™ d (i, 1', to) + t pd (i, V, m)5t™ d (i, I, m) 



(14) 



where t pd (i, I, m) represents the d-p transfer integrals be- 
tween I orbital of Co (l = d xy , d yz , d zx ) and m orbital of 
O (m = PxiPyiPz) at i site (i = 1 ~ 6), and 5tp d (i,l,m) 
is the change of t p d(i,l,m) due to a-mode phonon. In 
the absence of trigonal distortion, t p d(i, l,m) is given by 
Slater-Koster parameter t pd7r only, as shown in Fig. [3] 
For example, t pd {i = 1,1 = d xy ,m - p y ) = —t v d-K and 
tpd{i = 4, 1 = d xy ,m = p y ) = t p dTT. According to Harri- 
son's law [32]], t p dTr is proportional to a -4 , where a is the 
distance between Co and O. Then, 5t pd7T = —tpd^^Sa. 

Now, we calculate (5V BR )u. Since 5a = u BR /\/3 for 
breathing phonon, 5tpf'(i,l,m) = —tpd{i,l,m)^^u BR . 
For each I, t p d(i,l,m) is finite only for four sets of 
(i,m). After taking the summation of i and m in Eq. 
(dll), we obtain (5V$)u = -8t 2 



By transforming Eqs. 
we obtain 



5VT = 



I8|)-(p]l into (ai g , 



, e'g )-basis, 



f t SH1 



SH2 



<5U 



V2 




Ap d a V 2Mw SH 



(20) 

(21) 
(22) 



,BR 



pdir 

— 32t pd7T u BR / (\/3Apda) . By performing similar analysis, 
we obtain (5V BR ) lt > = tt 2 pd „u BR / {V3A pd a) for I ^ V by 
using the Slater Koster formula [H[, without necessity 
to use Harrison's law. As a result, the matrix represen- 
tation for breathing phonon in the (d xy , d yz , d zx )-basis is 
expressed as follows: 



Thus, the EPI originating from d-p transfer integrals is 
x 4 = represented by four coupling constants aj , aj, bj and 



SV 



BR 



3 -1 -1 
-1 8 -1 
-1 -1 8 



At 2 , 
p -u BR 



V3A 



pdCL 



(15) 



Equation (Ti~5"l) is transformed in the (ai g , e' 1 , e' 2 )-basis as 



8Vn 



• T 

; a 2 



BR 



aj 
aj 


f 2 



-BR 



24 



V3 A prf a V 2Afw B R 



(16) 



(17) 



In the same way, we can derive the EPI for shear 
phonons. For SHI mode (u SH i = ± 2t j=-(l, — 1, 0)), 

<5a = u/y/2 for Oi and O4, 5a = —uj\f2 for O2 and O5, 
and 5a = for O3 and 0$. Therefore, 8t pd (i, l,m) = 
-tpd{i,l,m)^u SYil for i = 1,4, t p(i (i, m)^-w SH1 for 
z = 2,5, and for i = 3, 6. For SH2 mode (ush~2 = 
±^(1,1,-2)), 5t^(i,l,m) = -t pd {i,l,m)^r a u SB:i 

for i = 1,2,4,5, and t pd (i, I, TO )^I uSH1 for i = 3,6. As 
a result, we obtain 



- SHI 

5V T = 



1 -1 

1 8 

-1 -8 



y/2t 



pdir 



,SH 



A 



pd ' 



5V 



SH2 



-16 

-1 8 
-1 2 



1 -1 



2 tpd-^ 

3 A pd a 



Sli 



(18) 



(19) 



B. EPI in the presence of trigonal distortion 



In mi A[ wc have calculated the EPI in the absence 
of trigonal distortion. However, C0O2 layer in Na^OoC^ 
becomes thinner due to the trigonal distortion, and it is 
increased by the water intercalation. Since this change 
of crystal structure can modify the EPI prominently, we 
calculate the EPI in the presence of trigonal distortion 
in this subsection. 

Owing to the trigonal distortion, the position of 
O ions move along the c crystal axis. Then, the 
changes of the coordinates of Oi, O2 and O3 are 
(— &/vo, — 6/vo, — 6/V3) and that for O4, O5 and Oq 
are (6/v3~, b/y/3, b/y/3), where b (> 0) is the magnitude 
of displacement. According to the neutron scattering 
measurement [39[, O1-C0-O5 angle 6*15 for unhydrated 
Na^CoC^ and hydrated Na^Co02 • J/H2O are 84° and 
82°, respectively. Thus, we can determine the value of b 
from #15 by solving the following equation. 



cos6>i E 



-(3m 2 - 2m)/(3m 2 - 2m + 1), 



(23) 



where m = b/(y3a). Using the values of b obtained in 
Eq. (|23[) . we derive the EPI in the presence of trigonal 
distortion as follows. 

First, we calculate the EPI originating from the change 
of the Coulomb potential. In the case of in — 0, the u a - 
linear terms of the Coulomb potentials are given in Eqs. 
([I])-©. Up to the first order of to = o/(V3a), they are 
modified as 
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^c BR (r) 



SVS m (r) 



^ c SH2 (r) 



2e 2 / 2V3 „ w A RR 

— + — 1 + 7m)(xy + yz + zx) u , 

47reo \ a z a* / 

2e 2 
47re 

2e 2 



9^2 
2^ 

3^ 



3\/2 



(yz — zx) J (1 + 7m)u 
V6 



SHI 



47re I 2a 4 

As a result, the coupling constants of, a!? , of and 02 are given by 

r 2e 2 2\/3 



(aT + y - 2z A ) + —r-{2xy - yz - zx) (1 + 7m)u' 



SH2 



47T£o a 2 V 2Mwbr 
2e 2 2y/3 



6? = 4^26. 



47T£o a 2 V 2Mo>br 
c _ 2e 2 2V3 



(l + 7m)), 
1 + ^(1 + 7m)), 
4 1(1 + 7 TO ). 



4tt£ a 2 V 2Mw S h 7a 2 
I 



(24) 
(25) 
(26) 

(27) 
(28) 
(29) 



The values of m at 9 15 = 84° and 82° are m ~ 0.051 
and 0.068, respectively. Therefore, of and 62 at 9 15 = 
84° and 82° are about 40 % larger than the values at 
#15 = 90°. The estimated values of of, a§, of and b% 
are shown in Table [IJ In calculating these values, we use 

= 14.4 cV (a = 2.0 A), M = 26.561 x 10~ 27 kg 
(the mass of 16 ion), r c i = 0.6 A(radius of Co ion [4C||). 
w B r = 580 cm -1 and w SH = 480 cm -1 [Hj]. In Eqs. 
(f24"|) - (|29j) . we show the results up to the first order of m. 
In calculating the values of a\, a-i, b\ and 62 in Tabic Q] 
we use the exact expressions for them with respect to to. 

Next, we calculate the EPI due to the change of d- 
p transfer integrals. In the absence of trigonal distor- 
tion, we have only to consider the contribution of t pc j w to 
t p d(i, I, m). In contrast, we have also to consider the con- 
tribution of tpda- in the presence of trigonal distortion, 
and the d-p transfer integral t n d( i, I, m) is given by us- 
ing the Slater-Koster formula [38|]. Then, we can obtain 
St" d (i, I, to) with the aid of the Harrison's law t cx a~ 4 
[37f . As a result, EPI can be derived from Eq. lfl4]). The 
obtained values of , aj, bj and oj are shown in Table 
HI In calculating these values, we use t v d a = —1-70 eV, 
t-gdj! = 0.785 eV and A p d = 1.81 eV according to Ref. 



' 1/dTV 

19]. 



015 


90° 


84° 


82° 




015 


90° 


84° 


82° 


a? 


-1.01 


-0.96 


-0.94 




T 
Ol 


-0.0675 


-0.121 


-0.147 


a? 


-1.09 


-1.08 


-1.07 




aj 


-0.101 


-0.162 


-0.189 


b? 


0.0621 


0.0874 


0.0971 




bl 


0.0557 


0.108 


0.136 


b'i 


0.0110 


0.0142 


0.0146 




bl 


0.0262 


0.0524 


0.0629 



TABLE I: Obtained electron-phonon coupling constants orig- 
inating from the change of the Coulomb potential (left panel) 
and the change of d-p transfer integrals (right panel) for 
0i5 = 90°, 84° and 82° 



015 


90° 


84° 


82° 


ai 


-0.215 


-0.218 


-0.218 


a.2 


-0.239 


-0.249 


-0.252 


bi 


0.118 


0.196 


0.233 


62 


0.0372 


0.0665 


0.0775 



TABLE II: Total electron-phonon coupling constants for 
015 = 90°, 84° and 82° 



electrons and optical phonons are represented as follows. 



We have calculated the EPI between ti g electrons and 
zone center phonons (q=0). Hereafter, we neglect the q- 
dcpcndcnces of the EPI matrix elements for simplicity, by 
expecting that their q-dependences are smeared out after 
the q-summation. The dimensionless displacement due 
to a-mode phonon in EPI can be expressed as the sum 
of creation and annihilation operators for phonon, = 
b a .q + b' a . Then, the Hamiltonian of EPI between t^g 



BR,SH1,SH2 

^pi = -7= £ E 4 +q .^ k(T (& Q , q + &L-J30) 

a k,q,cr 

where Cko- = (ck,oi g ,<7 c k.e^.o- c k.e«.o-) is the column vec- 
tor of annihilation operators for electrons, and o Q q is 
annihilation operator of a-modc phonon. Then, V a has 
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the following form: 



BR 



V 



ySKl = 



(11 











ft 2 











(12 





h 




h 





i) 





-b 2 





(31) 



where 



and bj 



ySH2_ 



b? 




The 



(i = l,2) 

estimated values of these coupling constants are shown 
in Table HT1 In estimating a\ and a-z, we have considered 
the screening effect by taking account of the chemical 
potential shift dp: To conserve the electron number, dp 
should satisfy the relation pd(0)(Sed + 6p) + p p (0)Sp = 0, 
where Pd(0) and p p (0) are the DOS of d-electrons and p- 
electrons at the Fermi level [28| . Then, the shift in the 3d 
level is 5id = fed + Sp = p p (0)/[pd{0) + pJOMSed- Since 
p p (0)/[p d (0) + p p {0)} « 0.2 in Na x Co0 2 [MOlI, both a x 
and 02 are reduced to 20% of their original (unscreened) 
values. On the other hand, such a screening effect is 
absent for shear phonons since Ti{V SH } is equal to 0. In 
later sections, we use b\ = 0.20 and b 2 — 0.067 (eV) for 
the unhydrated Na a; Co02 {0 15 = 84°), and b x = 0.23 and 
b 2 = 0.078 (eV) for the hydrated Na K Co0 2 (0i 5 = 82°). 
In both cases, we set a\ = 0.22 and 02 = 0.25 (eV). 



Both p a ig{uj) and p eg >(uj) in Eqs. (|33|) and (|34j) are equal 
to the values in the tight binding model at uj = and 
uj = A. Note that p alg (uj)du) = p eg > (u>)duj = 
0.8 in accord with the band calculation [lj] as well as 
the tight binding model [l!|. Here, the number of hole is 
0.65, which corresponds the hole number of Na .35CoO2. 



(a) 



(b) 






3- 


Peg' 






2- 




1 - 








La, 



A 



-1 



co(eV) 



III. DERIVATION OF GAP EQUATION 

In this section, we formulate the model and derive the 
linearized gap equation. First, we explain the energy de- 
pendence of the DOS in the present model. In Naa;Co02, 
there are three electronic bands (one ai ff -band and two 
e^-bands) near the Fermi level, which are composed of 3d 
orbitals of Co ion and 2p orbitals of O ion in C0O2 layer. 
In these bands, the ai g -band crosses the Fermi level, and 
the e^-bands are below the Fermi level (valence bands). 
Figure |4] shows the d-electron DOS of each band per spin 
that is obtained by the tight-binding model given in Rcf. 
[l9l |. In the numerical study, we approximate (i-electron 
DOS by the following rectangle function. 

! 1.38 ( -0.92 < uj < -0.73 ) 
0.38 ( -0.73 < uj < 0.3 ) , (33) 
3.65 ( 0.3 <oj< 0.34 ) 

{0.15 ( A -1.25 <uj< A -0.85 ) 
1.75 ( A -0.85 <uj< A -0.65 pf) 
0.6 ( A -0.65 <lj< A ) 

where A (< 0) is the top of the e' g bands measured from 
the Fermi level. The total d-electron DOS is p a i g (oj) + 
2p eg i(w). Here, we regard A as variable. Equations (|3"3"1) 
and (|34j) are shown in Fig. 5] (a) and (b), respectively. 



FIG. 4: (Color online) d-electron DOS for the ai s -band and 
the eg-band per orbit and spin. Solid line represents the 
DOS used in the present study (solid line), and dashed line 
represents the DOS given by the tight-binding model. Note 
that p eg i (uj) in the tight-binding model takes small values for 
ui > A, since the ai 9 -band that forms the large FS around T 
point contains small amount of e 9 -orbitals. 



Although LDA band calculation predicts A > (i.e., 
presence of e' g hole pockets) [ 141 , t he relation A < had 
been confirmed experimentally |15l . UtI [l8j . To solve this 
discrepancy, many people had studied the strong corre- 
lation effect, which is neglected in the LDA, using the 
Gutzwiller approximation |42ll and the dynamical mean 
field approximation [43|, |44|. References [13, [44| claim 
that e' g hole pockets predicted by the LDA disappears 
(i.e., A < 0) due to strong Coulomb interaction. We had 
also pointed out in Ref. [28[ that the EPI (due to shear 
phonon) shifts the e^-bands downwards by -0.05 eV. In 
the present paper, we treat A as a model parameter. 



A. gap equation 

In Ref. [28[, we have studied the linearized gap equa- 
tion for s-wave superconductivity in the 11-band d-p 
model. The obtained gap function (ie n ) does not have 
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the off-diagonal components (I ^ I') where I, I' are the in- 
dexes of orbital (I, I' = a lff , e'g 1 , e'g). Moreover, gap func- 
tions for e'g 1 and e'g are the same. Then, the linearized 
gap equation given by Eq. (8) in Ref. [28| , which is 3 x 3 
matrix equation, can be reduced to the following 1x2 
matrix equation: 



A0(ie„) 
Fi(ie n ) 



T^V cS (iui m )F(ie n - iu> m ), 

III 

^|G n ,(k,i £ „)| 2 0;(i^) 



(35) 



k./' 



= m 



(k n )l^|G„(k,k„)| : 



(36) 



where 0(ie„) = (0 o i fl (fc n ),0 efl '(ie n )) and F(ie n ) = 
(F a i g (ie n ), F eg i (i£„)) are column vectors of the gap func- 
tions and the linearized anomalous Green functions, re- 
spectively, and G;p(k, ie n ) is the normal Green function. 
s n = (2n + 1)ttT and u m = 2mirT are the Matsubara fre- 
quencies for fermion and boson, respectively. In Eq. (|36[) . 
the eigenvalue of the gap equation, A, increases mono- 
tonically as the temperature decreases, and T c is deter- 
mined by the condition A = 1. In Eq. (|36[) . we ignore 
the contribution oil ^ V since the off-diagonal compo- 
nents of Green function arc small. Then, the diagonal 
component of Green function is approximately given by 
G„(k,k„) = (k„ -El- Z l (k,ie n ))-\ where E'(k,i£n) 
is the d-electron self-energy. For the effective interaction 
of Cooper pairs V e e(iu) m ), we consider the EPI given by 
Eq. (|3T)|) and the Coulomb interaction Hu + Hj, where 



H, 



U_ 

N 



J2 H c k^T c k',ia Ck '-q' ; .i Ck +q.^T ( 37 ) 



k,k'.q I 



N 

k,k',q;,Z'^; 

where U and J are the on-site Coulomb interaction in 
the same orbital and the pair hopping potential between 
different orbits, respectively. Hereafter, we set J = U/10 
according to the first principle calculation in Ref. 
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Since we have assumed that p-electrons are noninteract- 
ing, V e ff(itu m ) works only for d-electrons. Note that we 
do not have to consider the Coulomb interaction between 
different orbitals (U 1 ) and the Hund's coupling term (Jh) 
since the off-diagonal components of gap function are ab- 
sent. Now, we treat Hjj and Hj in the Hartree-Fock 
approximation, by extending the theory of Morel and 
Anderson [45|. Then, V s{iuj m ) is given by 




26? A , 
b\ 2b\ ) DsH{lU>m) 



(39) 



where D a {iuj m ) = 2uj a /(ujl 



■ Li), 



,) is the phonon Green 



function for a-mode. ui a is the Debye frequency of a- 
mode phonon. The first two terms in Eq. (|39|) represent 
attractive forces due to EPI, and the last term represents 
the repulsive force due to the Coulomb interaction. Wc 
see that attractive force induced by breathing phonon, 
which corresponds to the first term in Eq. (|39[) . is re- 
duced by U (or U + J). On the other hand, attractive 
force due to shear phonon is reduced by J (or 2J). 

Here, we calculate Fi(\e n ) in Eq. (|36[) . Instead of 
summing k in the right-hand side, we perform the energy 
integration as 



1 fC 

-£|G„(k,k„)| 2 = / 



du, (40) 



where pi(u>) is the d-electron DOS defined in Eqs. 
(f3"3")) - (j3~4"|) , and Zi(w) is the renormalization factor 
given by the self-energy. The renormalization fac- 
tor due to EPI at w = is given by z/(0) = 



1 



;=0 



j=0 



, where Yif R (ui) and 

S^ H (w) are the self-energies given by breathing and shear 
phonons, respectively Since the renormalization due to 
the self-energy takes place only for \uj\ < tj a [28|, we 
approximate zi(u>) as 



I - 



LJ = 
BR/ 



(|cj| < ujsh) 

(w S H < |w| < W B r) 
("BR < M) 



(41) 



Note that T c is reduced by Zi(ui). In the one-loop ap- proximation [28|, T,f R (uj) and £f H (cj) are given by 

S™( ie ") = T^a 2 G a i s (ie„-^ m ) J D B R(^ m ),(42) 
m 

Se S R (*£„) = Tj24Geg(ie n -lLJ m )D BR {iLJ m ), (43) 

m 
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T^2b\G eg {ie n - iu m )D SIl (ioj m ),(AA) 
m 

{b\G a ig{ie n - iu m ) 

m 

+2b\G eg (ie n - iu> m )} D S}i (iuj m ). (45) 



In deriving above equations, we have put D a (iuj m ) = 
■£-9(tu a — \uj m \) for simplicity. Then, the differential coef- 
ficient of self-energy at u> — is expressed in the following 
form. 





Peg' {£) , 

Palg(s) 
Peg'(£) 



CO 



-de, 



BR 



-de. 



(46) 
(47) 



^SH 



In the next section, we solve Eqs. (p5"j) - (j56")l self- 
consistently using Eqs. (|4"0")) - (|41[) . and obtain T c that 
is determined by the condition A = 1 in Eq. (1351) . 



IV. NUMERICAL RESULT 

Here, we perform numerical study: In calculating T c , 
important model parameters are a\, b\ and A: a\ (bi) 
represents the EPI due to breathing mode (shear mode) 
in the ai 9 -orbital (between a\ g - and e^-orbitals), and A is 
the energy of the top of the e' g bands measured from the 
Fermi level in Fig. 0](b). Note that A < experimentally 
E3, El]- The values of ai and &i obtained in the 
present paper is slightly larger than those used in Ref. 
|28|. Therefore, the mass enhancement factors in the 
present model given by Eqs. (j46|) and ([47|) are slightly 
larger than those in Fig. 2 of Ref. [28|: In the present 
model, ^^(O) is around 2 for A ~ —0.02 eV, and it 
decreases moderately as A decreases. 

Then, T c is obtained by solving the set of gap equa- 
tions (|55"l) - ({56")) . In the absence of Coulomb interaction, 
the obtained T c -A phase diagram is qualitatively simi- 
lar to Fig. 4 in Ref. [28[, where obtained T c is about 10 
times higher than the experimental T c . In this section, we 
solve the gap equation in the presence of strong Coulomb 
interaction between d-electrons. We will show that the 
obtained T c is reduced to be comparable with experimen- 
tal T c by assuming reasonable values of U and J. In the 
presence of shear phonons (&i 7^ 0), gap function for e' 



orbital 



'eg' 



[ie n ) is finite even if the e' bands sink below 



the Fermi level (A < 0). Due to this interband transi- 
tion of Cooper pairs, which we call the valence band SK 
effect, s-wave T c is considerably enhanced. 

In Fig. [SJ we show [/-dependence of T c at A = —0.03 
eV for unhydrate (61 = 0.2 eV and b x = 0.067 eV) and 
hydrate (61 = 0.23 eV and b 2 = 0.78 eV) systems, re- 
spectively. We see that T c decrease monotonically with 
U since the attractive force due to phonons is reduced by 
U. In the absence of shear phonons (61 = b 2 = 0), T c 



is independent of A. Experimental T c ~ 4 x 10 -4 eV is 
realized at U ~ 1.75 eV, whereas T c is almost zero for 
U > 4 eV. In the presence of both breathing and shear 
phonons, obtained T c is considerably enlarged owing to 
the valence band SK effect induced by shear phonons, 
which is represented by b\ in Eq. For U ~ 4 eV, 

obtained T c for hydrate system is more than three times 
larger than that for unhydrate one. In §V A[ we explain 
the reason why SK effect is strong against large Coulomb 
interaction. 

A=-0.03 (eV) 



e 15 =82°(with water) 
e 15 =84°(without water) 
^ =b 2 =0 




[/(eV) 

FIG. 5: (Color online) [/-dependence of T c for A = 0.03 eV. 



In Fig. [6l we show the A-dependence of T c at U = 4 
and 6 eV for unhydrate and hydrate systems, respec- 
tively. Note that the cluster calculation using quantum 
chemical ab-initio methods suggests U = 4.1 ~ 4.8 eV 
[34| . In both cases, T c is almost zero if wc drop shear 
phonons (bi = 62 = 0). In the presence of shear phonons, 
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on the other hand, obtained T c is comparable to or higher 
than experimental T c due to the valence band SK ef- 



U=4 eV 



feet for |A| < 



In Fig. [6l T c increases mono- 



tonically with increasing A (< 0) due to SK effect, in 
particular for | A| < wsh- ARPES study in hydrated 

Na x Co0 2 -2/H 2 [H suggests that A 0.03 eV, which 

is closer to the Fermi level than the Debye frequency of 
shear phonon o>sh—480 cm" 1 ^0.06 cV. On the other 
hand, ARPES studies in unhydrated Na x Co0 2 [H, GJ] 
suggest A ~ —0.1 eV. Therefore, our results are consis- 
tent with these experimental results. 

The theoretical value of T c in Fig. [6] for |A| ~ 
0.02 — 0.03 eV is much larger than experimental T c , even 
for parameters of "without water". Therefore, we also 
calculate T c for smaller values of b\ and 6 2 in Fig. [6l 
&i = 0.18 eV and b 2 = 0.062 eV. The obtained T c for 
|A| ~ 0.02 — 0.03 eV is comparable with experimental 
value T c = 5K as 0.04 eV. 

In Appendix, we analytically study the transition tem- 
perature by taking account of the valence band SK effect. 
The obtained expression for T c is given in Eq. (IA14I) . 
where the effective coupling constant A e ff is given in Eq. 
(|A15|) . Since A 2 , A3 cx b\, the second term in Eq. (|A15|) 
represents the contribution due to the valence band SK 
effect. This term increases as A approaches zero since 
Aa is a decrease function of A (< 0). 

Finally, we summarize the results of this section. Ow- 
ing to the valence band SK effect, the obtained s-wave T c 
is relatively high even if we take account of the realistic 
Coulomb repulsion U = 4 ~ 6 eV. Moreover, not only A 
approaches zero but also EPI for shear phonons (61 and 
62) increase due to the change of crystal structure by 
water intercalation, as explained in ^11 BI For this rea- 
son, s-wave superconductivity is realized against strong 
Coulomb interaction in hydrated Naa,Co0 2 • j/H 2 0. 



DISCUSSION 



Why SK mechanism can overcome strong 
Coulomb interaction? 



In mVl we have shown that T c becomes at rela- 
tively small t/(~ 4 eV) when we consider only breath- 
ing phonon, while large T c is realized even for U > 4 
eV when we consider shear phonons as well as breathing 
phonon. This result indicates that the SK mechanism is 
suppressed by the Coulomb interactions only slightly. In 
fact, the attractive force due to shear phonons is reduced 
by J ~ U/10, which means that the reduction of the 
SK mechanism due to Coulomb interaction is very small. 
Moreover, J is prominently reduced by the retardation 
effect, as we will explain below: 

In a single band model, the retardation effect was dis- 
cussed by Morel and Anderson [45[ . In Appendix, we ex- 
tend their theory to the multi-band model for Na^CoO^ 
The renormalized Coulomb interactions are obtained by 
solving set of equations (|A7[) - (|A10[) . Here, we approxi- 
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r c (eV) 



0.005 



9 15 =82°(with water) 
9 15 =84 "(without water) 
^=0.18, £> 2 =0.062 
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FIG. 6: (Color online) A-dependence of T c at [7=4 eV (top) 
and !7=6 eV (bottom). The dashed and dotted line represent 
T c for hydrated Na a; Co02 -j/H 2 (61 = 0.223, b 2 = 0.0775) and 
unhydrated Na^OoO^ (&i = 0.196, 62 = 0.0665), respectively. 



mate the DOS for the a\ g - and e^-bands as p a i g (uj) = 
Pai g 0(W- \cj\) and p eg >(u) = p eg >0(A - oj)6(W + oj) [see 
Fig. [9] in appendix]. Here, 2W is the bandwidth of the 
ai g -band. In a realistic parameter regime U J, these 



12 



equations can be solved easily. The obtained results arc 



U, 



alg 



u; 



alg, eg' 



J. 



eg', eg 1 



U/(1 + AU), (48) 

U/(1 + BU), (49) 

J/(l + AU)(1 + BU), (50) 

J/(1 + SC/) 2 , (51) 



where U* lg (U* g , ) is rcnormalized Coulomb interaction 
between a lg electrons (e' g electrons), J* lg ^ eg , (Jt g ', eg ') 
is renormalized pair hopping between a\ g and e' g or- 
bitals {e' g orbitals). A = p a i g ln(W/tOD) and B = 
\p e g' Ixi(W/(Jd) are the renormalization factors. In the 
present model, A w 1.2 and B w 1.0 H|. (Note that 
Pai g U* lg {p eg 'U* gl ) corresponds to Anderson-Morel pseu- 
dopotential \x for a single band model.) Since the renor- 
malization factor for J in Eqs. (|50p and (|5ip is square of 
that for U, J* becomes very small. In the present model, 
U* lg /U « 1/7 and J^/J « 1/42 at U « 5 eV. For 
this reason, the effect of the Coulomb interactions on the 
SK mechanism becomes remarkably small. Therefore, in 
the case of |A| <C wd, s-wave superconducting state is 
realized owing to the SK effect even in the presence of 
strong Coulomb interaction. 



B. Isotope effect on T c 

Recently, Yokoi et al. studied the isotope effect on T c 
in Naa;Co02 • 2/H2O by substituting 16 O atoms with 18 O 
atoms [47J ■ They found that the isotope effect coefficient 
a (T c oc M~ a , where M is the mass of O ion) is consid- 
erably smaller than the simple BCS value 0.5. Here, we 
study the isotope effect on T c in the present model, and 
explain that the value of a approaches zero because of 
the strong Coulomb interaction. 

According to the BCS theory, T c in a single band 
model is T c = 1.13u D e-^ x " and A* = 2N(0)a 2 /cj D - 
where N(0) is the density of states, a is the EPI, and 
H* = N(0)U/(1 + N(0)U\u{W/lu d ) is the Anderson- 
Morel pseudo potential. In the case of fx* = 
(U = 0), 2N(0)a 2 /tu D is independent of M since a oc 
{Mujd)- 1 ' 2 oc M" 1 / 4 . Then, T c oc w D oc M' 1 ' 2 (i.e., 
a = 0.5). However, a becomes smaller than 0.5 in the 
case of U > 0, since /1* decreases as M (oc w^ 2 ) increases. 
In fact, a ~ in Ru and Zr, both of which are 4c?-clcctron 
s-wave superconductors. Since both N(0)U and lud/W 
are relatively large in Na^CoC^ • 2/H2O, the value of a is 
expected to be much smaller than 0.5. 

Figure [7] show the obtained isotope effect on T c at 
A = -0.025 cV and A = -0.04 eV, by assuming 
uj a oc M" 1 / 2 and 01,02,^1,62 oc Af -1 / 4 . For b\ and 
02 j we use the smaller values than the estimated ones in 
Table U to reproduce the experimental T c (~ 4.5 K). We 
use 61 = 0.18 eV, b 2 = 0.062 cV for A = -0.025 eV, 
and 61 = 0.20 eV, b 2 = 0.067 cV for A = -0.04 eV. 
(In both cases, we put ai = 0.22 eV and 02 = 0.25 cV, 
which are same with the estimated values in Table HIl ) 



For both parameters, the reduction of T c due to isotope 
effect, AT C , decreases with U. We see that AT C becomes 
almost at U ~ 5 eV, and it becomes positive at U = 6 
eV. Therefore, AT C in Naa;Co02 may be too small to 
probe it experimentally. At U — 3 eV, the obtained a 
is about half of the BCS value (a = 0.5). a becomes 
almost at U = 5 cV, and it becomes negative at U = 6 
eV (inverse isotope effect). Therefore, a approaches zero 
in Naa;Co02 -yH2 due to the strong Coulomb interaction, 
even if the s-wave superconductivity is caused by EPI. 
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FIG. 7: (Color online) Isotope effect on T c for at A = -0.025 
eV, 61 = 0.18 eV, b 2 = 0.062 eV (top), and A = -0.04 eV, 
61 = 0.20 eV, b 2 = 0.067 eV (bottom). The lateral axis shows 
the rate of ls O ions and the longitudinal axis shows T c /T c o 
where T c0 is T c for ie O ion. Note that 01 = —d\aT c /d\nM is 
the isotope effect coefficient. 
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C. Origin of the anisotropy of the superconducting 

gap 

In this paper, we have studied the s-wave supercon- 
ducting state due to EPI in Na^OoO^ • 2/H 2 0. However, 
the nuclear relaxation ratio 1/T± is proportional to T 3 
below T c , which suggests that the nodal superconduct- 
ing state is realized Presence of gap anisotropy is 
also indicated by the specific heat measurements below 
T c [H, 51] • In the present calculation, we have neglected 
the k-dependence of the superconducting gap for sim- 
plicity. Here, we discuss the possibility of realizing an 
anisotropic s-wave state in Na^CoC^ • yH 2 0, resulting 
from the coexistence of the strong EPI and the AF fluc- 
tuations. 

In (Y,Lu)Ni 2 B 2 C, strongly anisotropic superconduct- 
ing state is realized below T c ~ 15 K [50, [H|. The super- 
conducting gap becomes isotropic by introducing small 
amount of impurities, whereas T c is almost unchanged 
[52| . They are hallmarks of anisotropic s-wave supercon- 
ducting state. According to Ref. [50], the ratio of gap 
anisotropy reaches A max /A m i n ~ 10 in a clean sample. 
Recently, one of the authors of this paper had studied the 
mechanism of this anisotropic s-wave superconductivity 
[53|: By solving the strong coupling Eliashberg equation, 
he found that the s-wave superconducting gap due to 
the EPI can be strongly anisotropic even in a single FS 
model, if strong AF fluctuations exist. In this case, pairs 
of gap minima appear at points on the FS which are con- 
nected by the nesting vector Q. In the normal state of 
Na^CoC^ • 2/H2O, prominent AF fluctuations had been 
observed by 1 jT\ measurements [HI, [54| . According to 
Ref. [13], Na^OoC^ • yH 2 locates in the close vicinity of 
the AF quantum critical point. Recent neutron diffrac- 
tion measurement reports that the wavevector of the AF 
fluctuations is Q 2 = fM = (0,2tt/\/3) and (tt, ±7r/\/3) 

According to the study using the FLEX approximation 
[l9| . there are two candidates for the nesting vectors in 
Na x Co0 2 for A < 0. One is Qi = FK = (4tt/3,0) and 
(±27r/3, 27r/\/3) that originates from the a\ g large FS 
around T-point. The another one is Q 2 that originates 
from the nesting of e' g hole pockets near K-points. Note 
that the AF fluctuations with Q 2 appear when the top 
of the e' g bands is very close to Fermi level [l{|. When 
the hole number is 0.65, which corresponds to the hole 
number of Nao.35Co0 2 -1.3H 2 without oxonium ion, the 
resultant minimum gap points are shown in Fig [8] (a), 
when the wavenumber of AF fluctuations is either Qi or 
Q 2 . Note that the sign of the gap function is the same 
everywhere. On the other hand, when the hole number 
is about 0.5, which can be realized in the presence of 
oxonium ion as a substitute for water molecule (35[, the 
resultant minimum gap points are shown in Fig [8] (b) 
when the wavenumber of AF fluctuations is Q 2 . In both 
cases, anisotropic s-wave superconductivity is realized in 
the e^-band FS in Na x Co0 2 ■ yH 2 0. 



M 




FIG. 8: Expected anisotropic s-wave superconducting gap for 
(a) rih ~ 0.65 and (b) ~ 0.5, where rih is the hole number 
of t2 g orbital. The dashed lines represent the magnitude of 
gap function, whose sign is positive everywhere . Tops of the 
eg-bands, which are just below the Fermi level, are shown by 
dotted lines in (b). 



In Ref. [HH, we have shown that the relation l/T\ oc 
T 3 below T c can be realized in the anisotropic s-wave 
state under the influence of strong AF fluctuations. 
Moreover, we have recently verified that the two-gap type 
specific heat observed in Refs. [H, [49[ can be realized in 
this anisotropic s-wave state, by considering only the a\ g 
FS (56[. It is an important future problem to reproduce 
the anisotropic s-wave state microscopically based on the 
d-p Holstein Hubbard model for Na^CoO^ 



VI. SUMMARY 

In Na^CoC^, existence of the strong EPI is suggested 
by ARPES measurements [Hj], and the absence of im- 
purity effect on T c strongly indicates the realization of 
the s-wave superconducting state without sign change of 
the gap function. In the present paper, we have studied 
the electron-phonon mechanism of superconductivity by 
considering two relevant optical phonon modes (breath- 
ing and shear phonons), and found that strong pairing 
interaction is caused by the interband hopping of Cooper 
pairs induced by shear phonons. This mechanism is im- 
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portant even if the top of e' g electron band is close to 
but below the Fermi level as suggested experimentally. 
Therefore, this valence-band SK mechanism is the origin 
of s-wave superconductivity in Na^OoC^, overcoming the 
strong Coulomb interaction U ~ 5 eV. 

In this paper, we have derived the EPI for breathing 
and shear phonons by considering both the Coulomb po- 
tential and the transfer integrals. The estimated EPI for 
shear phonon is prominently increased by water interca- 
lation, resulting from the increase of trigonal distortion 
of C0O2 layer. According to the point charge model, the 
top of the e' g bands A(< 0) is expected to approach the 
Fermi level due to water intercalation. Both effects in- 
duced by the water intercalation will raise T c due to the 
SK mechanism. 

Based on the obtained model Hamiltonian, we deter- 
mine T c by solving the strong coupling Eliashberg equa- 
tion. The SK mechanism is seldom damaged by the 
Coulomb interaction since the pair hopping J, which 
is the depairing force for the SK mechanism, is much 
smaller than U. For this reason, experimental T c ~ 5 K 
is realized irrespective of the realistic Coulomb interac- 
tion U ~ 5 eV. Wc have also studied the oxygen isotope 
effect ( 16 -^ 18 O) on T c . In the absence of Coulomb 
interaction, T c decreases by the isotope substitution in 

— 1/2 

proportion to ujd oc Mq . However, we found that the 
isotope effect on T c becomes very small for U = 4 ~ 6 eV, 
since the renormalized Coulomb interaction (Anderson- 
Morel potential) is reduced with the decrease of 10 d- 
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APPENDIX A: ANALYTICAL EXPRESSION FOR 
THE TRANSITION TEMPERATURE IN 

Na x Co02 yH 2 

In this appendix, we derive the analytical expression 
for s-wave T c in Naa;Co02 • 1/H2O. For this purpose, we 
simplify the model further. Here, we approximate the 
DOS for the a\ g - and e^-bands as p a i g {uj) = p a igO(W ~ 
\uj\) and p eg '(uj) = p eg '9(A — uj)9(W + w), which are 
shown in Fig. Hereafter, we promise that A < 0. We 
also approximate the phonon Green function by the step 
function D(iu> m ) = ^■6'(ojd — |^m|), and assume that the 
Debye frequencies of breathing and shear phonon are the 
same for simplicity (wd = wbr = ^sh)- Then, the gap 
equation at T = T c given in Eq. ([55]) becomes 



(i£n) 



E 

|"m|<WD 



2b\ 



— F(ie„ - iw TO ) 



U 2.J 
J U + J 



F(i£ T 



(Al) 



Here, we assume that (f>(is n ) is given by the step function 



<^(i£n) 



0" (|e„| > wd) 



(A2) 



First, we consider the case of \e n \ *C ojb- Then, the gap 
equation for (j) 1 is obtained from Eq. (|A1[) . 




Y?\GeA^ t )\^\ g , 



(A3) 



and 



where and E" represent § EkE| £(! |< WD 

Ek E| £j! |>w D ! respectively. Next, we consider the case 
of |e„| > w D , where F(ie n - iuj m ) oc ((f>™ lg , <j% g , ) for 



< ujd. Then, the gap equation for 





E n |G al9 (k,i £ ,)| 2 ^ L 
£ n |G es ,(k,i e ,)' 2 ^ 



£ I |G efl ,(k,fc / )| 2 <' 



eg' 



(A4) 
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Here, wc ignored the phonon contribution that is propor- 
tional to 4> , since it is much smaller than the first term 
in Eq. jA4}. 




T* 

u alg,eg' 



J JBU: g , UAJ* alg>eg , 



-JBJ, 



eg .eg ' 



^ eg 1 , eg 1 ~ J a TBI) eg i J AJ a i geg , 

— UBJ eg , eg , , 



(A9) 
(A10) 



where A = \G a ig(k, iei)\ 2 and B = 

^2 |G es '(k, ie^)! 2 . If we assume |A| <C ujb, then A = 



Palg ln(W/o; D ) and B = Peg , Q I(| + lArctan(f ))de 

\p eg > \n{W/u>T>). 



We now calculate T c from Eq. (|A5|) . In the case of 

|A| <C c^d "C W, ^2 \Ge(k,ie n )\ 2 is approximated as 
follows. 



FIG. 9: (Color online) The simplified DOS for the ai g - and 



e 3 -bands. 



Using eq. (|A4[) . cj> u in Eq. (|A3[) can be eliminated. 
Then, the obtained gap equation for <p l is given by 



oil' 



efi 



4 



E'lGaigiKiee)] 2 ^ 



(A5) 



2b 2 



b\ a 



261 



U, 



alg 



j 

2 ^alo,e 



.'/■'</ 



alg,eg' ^ eg' 



J, 



(A6) 



eg', eg' 



where V* s is the renormalized effective interaction. U* lg 
(U* g ,) is the renormalized Coulomb repulsion for the a\ g 
(e' g ) orbital, and J* lg eg , (J* g , ,) is the renormalized 
pair hopping between (a\ g , e' g ) orbitals ((e^, e' g ) orbitals). 
These renormalized Coulomb interactions satisfy the fol- 
lowing simultaneous equations |46| . 
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As a result, according to Eq. (|A5|) . the equation for T c 
is obtained as 



A* In + g 



A* In 



A 



n + ln 143^ 



(A13) 



where A} = (g - U* lg ) z aigPaig , \* 2 = (g - 2J a ^ e ,) z e ,p e , , A3 = 

1 



2(q|+2o^ 



U e> - J e' e< Z e'Pe' and A A 
s g' g/ s s 



^ 2 111 V|A| 



7^7 — J a lg ,e' J Z ai s Pai s i A 4 — 



By solving eq. (|A13|) , T c is given by 



T c = 1.13ciJDCxp 



An Aq 



In Eq. (|A15[) . the first term originates from breathing 
phonon, and the second term originates from valence 



(A14) 
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band SK effect. Aa monotonically decreases with increas- 
ing A. Therefore, T c increase as the top of the e' g band 



approaches the Fermi level. 
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